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ABSTRACT 

We use the high angular resolution in the near-infrared of the WFC3 on HST to determine YHVz 
color-color selection criteria to identify and characterize 1.5 < z < 3.5 galaxies in the HUDF09 and 
ERS (GOODS-South) fields. The WFC3 NIR images reveal galaxies at these redshifts that were 
undetected in the rest-frame UV HUDF/GOODS images, as well as true centers and regular disks 
in galaxies classified as highly irregular in rest-frame UV light. Across the 1.5 < z < 2.15 redshift 
range, regular disks are unveiled in the WFC3 images of ~25% of both intermediate and high mass 
galaxies, i.e., above 10 10 M@. Meanwhile, galaxies maintaining diffuse and/or irregular morphologies 
in the rest-frame optical light — i.e., not yet dynamically settled — at these epochs are almost entirely 
restricted to masses below IO^Mq. In contrast at 2.25 < z < 3.5 these diffuse and/or irregular 
structures overwhelmingly dominate the morphological mix in both the intermediate and high mass 
regimes, while no regular disks, and only a small fraction (~25%) of smooth spheroids, are evident 
above 1O 11 M . Strikingly, by 1.5 < z < 2.25 roughly 2 out of every 3 galaxies at the highest masses 
are spheroids. In our small sample, the fraction of star- forming galaxies at these mass scales decreases 
concurrently from ~60% to ~5%. If confirmed, this indicates that z ~ 2 is the epoch of both the 
morphological transformation and quenching of star- format ion which assemble the first substantial 
population of massive ellipticals. 

Subject headings: Galaxies - galaxies: formation - galaxies: evolution - galaxies: high-redshift 



1. INTRODUCTION 

The redshift interval 1.5 < z < 3.5 represents an im- 
portant epoch for the study of galaxy formation and 
evolution, encompassing both the peak of cosmic star 
formation (Lilly et al. 1996; Madau et al. 1996, 1998; 
Hopkins & Beacom 2006; Bouwens et al. 2007; Tresse 
et al. 2007) and the peak of active galactic nuclei ac- 
tivity (Warren, Hewett, & Osmer 1994; Fan et al. 2001; 
Croom et al. 2004; Barger et al. 2005). One efficient 
method for identifying large samples of galaxies at these 
redshifts is via application of color-based selection cri- 
teria to optical-to-near-infrared (NIR) imaging surveys. 
Established selection criteria now exist for the identifi- 
cation of z ~ 3 systems via the Lyman break in ob- 
served UV-to-optical colors (cf. LBGs: Steidel et al. 1996, 
2003), and 1 < z < 3.5 systems via the Balmer break or 
4000A break in observed optical-to-NIR colors (cf. EROs: 
Thompson et al. 1999; Franx et al. 2003; McCarthy et 
al. 2004; Van Dokkum et al. 2006; BX/BMs: Steidel et 
al. 2004; Adelberger et al. 2004; and BzKs: Daddi et 
al. 2004, 2005). Importantly, with optical-to-NIR color- 
selection criteria it is possible to sample uniformly both 
actively star-forming and passively-evolving systems at 
high redshift (e.g. Daddi et al. 2004). 

The new Wide Field Camera 3 (WFC3) on the Hub- 
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ble Space Telescope (HST) has greatly enhanced the UV 
and NIR imaging capabilities of this observatory, offer- 
ing sharper imaging and greater throughput with a new 
complement of narrow- and broad-band filters. In partic- 
ular, the unprecedented depth and sharpness of the NIR 
WFC3 images further enhances our capacity to detect 
LBGs at the era of peak star formation in the universe — 
as well as at much higher redshifts; see, e.g., the WFC3 
detections of numerous galaxy candidates at z ~ 7-8 
(Bouwens et al. 2010a,b; Oesch etal. 2010; Wilkins et al. 
2010; McLure et al. 2010; Bunker et al. 2010; Finkelstein 
et al. 2010), and even z ~ 10 (Bouwens et al. 2011). 

In this paper we use the WFC3 data acquired so far 
in the context of the Hubble Ultra Deep Field 2009 
(HUDF09) and Early Release Science (ERS) programs 
to (i) extract i^i6o~band based galaxy catalogs down 
to, respectively, AB=27 and 25 mag; (n) derive optimal 
color-selection criteria ( YHVz) for the identification and 
characterization of galaxies at 1.5 < z < 3.5, based on 
the WFC3 NIR filter set (complemented by optical ACS 
data), and (ill) use the constructed WFC3 i^i6o-band 
source catalogs to investigate, through qualitative com- 
parisons with archival ACS optical data in the HUDF09 
and ERS fields, differences and similarities in the rest- 
frame UV and optical morphologies of 1.5 < z < 3.5 
galaxies. The efficiencies of the YHVz selection crite- 
ria presented here are established via comparisons with 
the spectroscopic and photometric redshifts of galaxies in 
both datasets. Alternative color-selection criteria, VIH 
and JHVI, relevant to the filter set of CANDELS/Wide, 
are also explored using the ERS dataset. 

The outline of the paper is as follows. In Section [2] 
we describe the imaging data now available from the 
HUDF09 and ERS observational programs, and briefly 
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review the other datasets from which additional UV-to- 
IR magnitudes and spectroscopic redshifts are obtained. 
In Section [3] we detail the construction of our master 
source catalogs based on the HUDF09 and ERS Hiqq 
imaging, and the matching of these catalogs to multi- 
wavelength photometry from previous studies. In Sec- 
tion [4] we describe the estimation of photometric red- 
shifts, specific star formation rates, and stellar masses 
via spectral energy distribution (SED) fitting. In Section 
[5] we present YEVz color-selection criteria for the iden- 
tification of galaxies at 1.5 < z < 3.5 using the WFC3 
and ACS filters. In addition, we carefully quantify selec- 
tion efficiencies and contamination rates for these crite- 
ria according to the spectroscopic, and (up to) 15 band 
photometric, redshifts of galaxies in our master source 
catalogs. Furthermore, the characterization of passively- 
evolving galaxies at 1.5 < z < 3.5 via a simple ex- 
tension of our s elect ion criteria (passive YHVz) is ex- 
plored in Section [5^4} and the alternative color-selections 
JHVI and VIE, using filters from the CANDELS/Wide 
program (excluding Y105 which is available only within 
the GOODS-N field; Koekemoer et al. 2011), are pre- 
sented in Section [53] In Section [6] we publish our master 
source catalogs containing positions, photometric aper- 
tures, and photometric redshifts for all robust detections 
to #i 6 o < 27 mag in the HUDF09 and H 160 < 25 mag 
in the ERS. In Section [7] we compare the rest-frame UV 
and optical morphologies of galaxies in the HUDF09 and 
ERS fields at 1.5 < z < 3.5, and examine the mass scales 
over which key structural sub-types are revealed by the 
WFC3. We summarize in Section [5J Note that we adopt 
a cosmological model with Cl\ = 0.75, Qm — 0.25, and 
h = 0.7, and that all magnitudes are quoted in the AB 
system throughout. 

2. DATA 

2.1. HUDF09 

The HUDF09 program (Bouwens et al. 2010) now in 
progress will observe a ~4.7 arcmin 2 central region of 
the HUDF (Beckwith et al. 2006) with the new Wide 
Field Camera 3 (WFC3) on the Hubble Space Telescope 
(HST) for a total of 96 orbit^J The exposure time is 
divided between three NIR filters: Y 105 (F105W), Ji 25 
(F125W), and H 160 (F160W). The first epoch data used 
here total 60 orbits (16 in Y105, 16 in J125, and 28 in 
#160)5 an d reach a uniform depth in all filters of ~28.8 
mag in a 0.2 arcsec radius aperture (despite the neces- 
sary rejection of two orbits of Y105 imaging due to severe 
persistence issues). Our reduction of the WFC3 imag- 
ing data using the multidrizzle package (Koekemoer et 
al. 2002) produces final science grade images with PSF 
FWHM of ~0.15 arcsec and a pixel size of 0.06 arcsec 
(see Bouwens et al. 2010). Rebinning of the HUDF ACS 
£435 (F435W), F 60 6 (F606W), i 775 (F775W), and z 850 
(F850LP) images to the scale of our WFC3 frames was 
also performed to facilitate aperture-matched photome- 
try. 

2.2. ERS 



6 A further 96 orbits in the HUDF09 program will be dedicated 
to imaging, at similar depth, of the two nearby HUDF05 fields 
(Oesch et al. 2007, 2009). 



The Early Release Science (ERS) program has ob- 
served a northern section of the GOODS-South field (Gi- 
avalisco et al. 2004) in the WFC3/IR channel in 10 point- 
ings for a total of 60 orbits (cf. Windhorst et al. 2010). 
The exposure time was divided evenly between three NIR 
filters: Y 098 (F098M), J 125 (F125W), and E 160 (F160W), 
reaching a depth of ~27.5 mag (in a 0.2 arcsec radius 
aperture) over an area of ~40 arcmin 2 . Our reduction 
of this data (Bouwens et al. 2010) was performed using 
the same procedures as for the HUDF09 WFC3 imaging, 
resulting in final science grade images with PSF FWHM 
~0.15 arcsec and a pixel size of 0.06 arcsec. The GOODS 
ACS £435 (F435W), V 606 (F606W), i 775 (F775W), and 
^850 (F850LP) images were again rebinned to the same 
scale as the reduced WFC3 frames. 

2.3. Ancillary Datasets 

A wealth of complementary data, including pho- 
tometry at wavelengths beyond the range of the 
HUDF/HUDF09 and GOODS/ERS HST imaging has 
been gathered for objects in the HUDF and GOODS- 
South fields by a number of different teams. Much of 
this data where publicly available has been compiled in 
the GOODS-MUSIC catalog (Grazian et al. 2006; San- 
tini et al. 2009), which we employ in this study for ob- 
jects with unambigu ous m atches in our primary source 
catalogs (see Section [D] below) . The GOODS-MUSIC 
catalog contains PSF- and aperture-matched photome- 
try (derived via the ConvPhot package; De Santis et al. 
2007) in the following additional filters: ESO 2.2-WFI 
and VLT-VIMOS U, ESO VLT/ISAAC J, JET, and K 8 , 
Spitzer IRAC 3.6, 4.5, 5.8, and 8.0 /mi, and MIPS 24 /mi. 
A compilation of spectroscopic redshifts is also included 
in the GOODS-MUSIC catalog, mainly sourced from the 
GOODS (Vanzella et al. 2005, 2006, 2008), K20 (Mignoli 
et al. 2005), and VVDS (Le Fevre et al. 2004) projects. 

In addition to the GOODS-MUSIC catalog, we also 
employ the published spectroscopic redshifts and spec- 
tral classifications for passive galaxies at z ~ 2 in the 
HUDF and GOODS-South fields from Daddi et al. (2005) 
and Cimatti et al. (2008) (hereafter pBzK). The for- 
mer sample consists of 7 high redshift sources in the 
HUDF with passive BzK colors, confirmed via analy- 
sis of low resolution GRAPES spectra, of which 3 lie 
in the region of WFC3/IR imaging from the HUDF09 
program. The latter sample consists of 13 galaxies from 
the GMASS catalog of Spitzer 4.5 /im detected sources 
in the GOODS-South field with a prominent Mguv fea- 
ture in their ESO VLT/FORS2 spectra, of which 10 lie 
in the region of WFC3/IR imaging from the ERS pro- 
gram. In total, 2 and 8 of these pBzK galaxies in the 
Daddi et al. and Cimatti et al. samples, respectively, lie 
within the target range of our YEVz color-selection cri- 
teria (1.5 < z < 3.5). A further 3 and 6 GMOS spectro- 
scopic redshifts from the 1 < z < 2 Roche et al. (2006) 
ERO sample were employed for objects in the HUDF09 
and ERS, respectively. Although only one of the corre- 
sponding objects was missing a 'confident' spectroscopic 
redshift in the GOODS-MUSIC catalog, the Roche et al. 
redshifts were preferred due to the higher confidence and 
precision of these measurements. Importantly, for one 
object (ERS0325) assigned z spec = 0.481 in the GOODS- 
MUSIC catalog, the value of z spec = 2.017 reported by 
the Roche et al. team is much more consistent with its 



observed UV-to-IR SED. Finally, for the purposes of star- 
galaxy separation we also refer to the PEARS-S star cat- 
alog of Pirzkal et al. (2009). 

3. EXTRACTION OF THE H 160 HUDF09 AND ERS SOURCE 
CATALOGS 

3.1. Details 

Master source catalogs based on all robust #i6o de- 
tections to #160 < 27 mag and #160 < 25 mag 
in the HUDF09 and ERS imaging, respectively, were 
constructed as follows. First, SExtractor (Bertin & 
Arnouts 1996) was run on each of the reduced #160 sci- 
ence images with four different choices of deblend param- 
eters representing hot (nthresh = 64), medium (nthresh 
= 16), cold (nthresh = 8), and very cold (nthresh 
= 2) extractions. In each case, the RMS frames out- 
put from multidrizzle were scaled to reflect the true 
signal-to-noise of the reduced science images before their 
application as weight maps, and the following key con- 
trol parameters were adopted: detect_thresh = 1.5cr, 
detectjninarea = 12, and min_cont = 0.0001. Based 
on these #160 extractions we then re-ran SExtractor 
in dual image mode to recover aperture-matched fluxes 
and Kron magnitudes from the remaining optical-to-NIR 
HST images (#435 to J125). Careful visual inspection of 
the output Kron ellipses of sources identified in each ex- 
traction was then performed to remove spurious detec- 
tions, and to identify the most appropriate aperture for 
flux measurement of each real source. When choosing the 
best aperture we considered both the positions of possible 
matches in the GOODS-MUSIC catalog and the appear- 
ance of each source in all seven ACS plus WFC3 (#435 
through #160 ) filters. In ambiguous cases we favoured 
the deblending of objects into distinct concentrations of 
Hiqq flux (the filter expected to best trace the underly- 
ing stellar mass distribution, modulo some uncertainty 
regarding the role of TP-AGB stars). For calibration of 
the NIR magnitudes we employed the latest (AB sys- 
tem) WFC3/IR zero-points available from the STSCi 
webpaga^J and for the optical magnitudes we adopted 
the ACDzero-points determined by the GOODS tearrR 
A total of 1052 and 3078 sources were thereby recovered 
to #i6o < 27 mag and H 160 < 25 mag in the HUDF09 
and ERS imaging, respectively. To define the master 
source catalogs employed in this analysis we limit our 
samples to those systems with high S/N coverage in all 
7 available ACS plus WFC3 filter^] and lying within the 
region of sky covered by the GOODS-MUSIC catalog, a 
total of 993 and 2630 objects in the HUDF09 and ERS, 
respectively. 

For all objects in our master source catalogs we 
searched for counterparts in the GOODS-MUSIC cat- 
alog (Santini et al. 2009) with matching positions and 
photometric apertures. All object pairs with centroid 
differences of <0.24 arcsec (i.e., 4 pixels in our drizzled 
WFC3 images) and #435, Vqoq, 2775, and zgso magnitudes 
differing by less than 0.5 mag between the two catalogs 
were classified automatically as successfully matched — 
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7 URL(http://www. stsci.edu/hst/wfc3/phot_zp_lbn) 

8 URL ( http : //archive . st sci . edu/pub/hlsp/goods/v2/ ) 

9 Identified as those objects occupying areas of the reduced sci- 
ence images with local background RMS less than twice the global 
background RMS. 
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Fig. 1. — The distributions of galaxies, spectroscopically- 
confirmed stars, and candidate stars in the if 160 apparent 
magnitude— Kron size plane in our master source catalogs to 
if 160 < 27 mag in the HUDF09 (top row) and if i 6 o < 25 mag 
in the ERS (bottom row). The red hashed area in each case de- 
notes the region within which small size outliers were classified as 
candidate stars unless resolved structure w as c learly evident upon 
visual inspection (as described in Section |3.2| ). The black, dark 
blue, and light blue (smoothed) contours respectively enclose 68%, 
95%, and 99% of the detected objects in each survey. 

426 sources in the HUDF09 and 1954 in the ERS fields. 
The remaining object pairs not satisfying these criteria 
were visually inspected, and a further 44 and 294 matches 
were thereby recovered for the HUDF09 and ERS fields, 
respectively. 

For most galaxies at #160 > 24.5 mag no counterpart 
could be identified due to the brighter GOODS-MUSIC 
selection limits (zs5o < 26 mag or K s < 23.5 mag). In 
total, 523 and 382 objects in our master source cata- 
logs from the HUDF09 and ERS, respectively, were ul- 
timately deemed to be missing valid counterparts in the 
GOODS-MUSIC catalog. 

3.2. Star- Galaxy Separation 

Star-galaxy separation in our HUDF09 and ERS sam- 
ples was performed as follows. First, all sources matched 
to counterparts in the GOODS-MUSIC catalog with 
spectroscopic redshift quality flag <2 and z spec = (2 ob- 
jects in the HUDF09, and 29 in the ERS) and all objects 
with centroids matched within 0.36 arcsec to members 
of the PEARS-S star catalog of Pirzkal et al. (2009) (an 
additional 3 and 16 sources in the HUDF09 and ERS, re- 
spectively) were identified as spectroscopically-confirmed 
stars. In addition, all ultra-compact outliers from the ob- 
served #160 apparent magnitude— Kron radius distribu- 
tions at log 10 (i?Kron/") < -0.125 x (H 160 - 20) + 7 (with 
7 = 0.2 for the HUDF09 and 0.075 for the ERS) not al- 
ready identified as spectroscopically-confirmed stars were 
visually-inspected and classified as candidate stars (2 in 
the HUDF09, and 35 in the ERS) unless resolved struc- 
ture was clearly evident. Figure [l] shows these candidate 
stars on the #160 apparent magnitude— Kron radius di- 
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Fig. 2. — Evaluation of the accuracy of our ZEBRA photometric redshifts for galaxies in the HUDF09 (top row) and ERS (bottom row). 
The panels in the left most columns present the one-to-one comparison between our photometric redshifts and available spectroscopic 
redshifts, while those in the colums second from the left evaluate our estimates of the photometric redshift uncertainties for these galaxies 
against the same spectroscopic benchmarks. (The green, yellow, and red shaded regions here mark deviations of < 5%, < 20%, and > 20% 
in Az/(1 + z), respectively.) The panels in the two right most columns illustrate the redshift distributions of galaxies in each field both as 
an ensemble, and after subdivision into those systems with spectros copi c redshifts, photometric redshifts with matched GOODS-MUSIC 
photometry, and HST-based photometric redshifts only. See Section [Xl] for further details. 



agnostic diagram. We exclude these sources from our 
analysis of galaxy color-selection efficiencies in Section [5j 

4. ESTIMATES FOR REDSHIFTS, SPECIFIC STAR 
FORMATION RATES, AND STELLAR MASSES 

4.1. Photometric Redshift Estimates 

We have computed photometric redshifts for all objects 
in our HUDF09 and ERS source catalogs via template 
fitting to the observed spectral energy distributions in 
the HST ACS and WFC3 filters (£435 to H 160 ), supple- 
mented by the additional photometry in the GOODS- 
MUSIC catalog (for a total of 15 passbands) where a 
valid counterpart could be identified. The actual tem- 
plate fitting was performed using the Zurich Extragalac- 
tic Bayesian Redshift Analyzer (ZEBRA; Feldmann et 
al. 2006). To recover optimal photometric redshifts ZE- 
BRA modifies an empirical template library (Polletta et 
al. 2007; plus model SEDs for very blue objects from 
Bruzual & Chariot 2003) to account for systematic mis- 
matches against the observed spectral energy distribu- 
tions. 

We quantify empirically our mean photometric red- 
shift uncertainties by comparison against matched ob- 
jects with confident (quality flag <2) spectroscopic red- 
shifts in the GOODS-MUSIC, Daddi et al., Cimatti et 
al., and Roche et al. catalogs (86 in the HUDF09, and 
446 in the ERS), as illustrated in the panels in the left 
most columns of Fig. [2] We thereby estimate typical un- 
certainties (Az = \z spec — Zphot |) °f ^ z ~ 0.063(1 + z) 
and Az ~ 0.039(1 + z) for the HUDF09 and ERS master 
source catalogs, respectively. The larger uncertainties in 
the HUDF09 photometric redshifts reflect the relatively 
small number of confident spectroscopic measurements 



available in this field for calibration of our SED template 
library. We also note low rates of catastrophic failures for 
the ZEBRA photometric redshift pipeline, namely that 
<5% of galaxies in the HUDF09, and <9% of galaxies 
in the ERS, are in disagreement with their spectroscopic 
redshifts by Az > 0.2(1 + z). 

At z < 1.5, where the greatest numbers of spectro- 
scopic redshifts are available (78 in the HUDF09, and 
355 in the ERS), the mean uncertainties are reduced 
to Az - 0.058(1 + z), and Az - 0.035(1 + z), respec- 
tively. Nevertheless, at z > 1.5 our photometric redshift 
accuracy remains satisfactory with uncertainties of only 
Az - 0.069(1 + 2), and Az - 0.056(1 + 2), respectively— 
although we note that in the case of the HUDF09 the reli- 
ability of this comparison is limited by the small number 
(8) of spectroscopic redshifts in this interval. Hence, as 
a further stage of quality control on our ZEBRA photo- 
metric redshifts we also evaluate here the performance 
of the output confidence intervals derived from the pos- 
terior probability distribution in redshift space for each 
galaxy. Specifically, we examine the ratio between Az 
(i.e., the absolute difference between the maximum like- 
lihood photometric redshift and the true spectroscopic 
redshift) and the width of the one sigma confidence in- 
terval for each galaxy (cf. Ilbert et al. 2009). If our poste- 
rior probability distributions are well estimated then we 
would expect the value of this quotient to be less than 
one for roughly 68% of each sample. In fact, our inves- 
tigation reveals a slight underestimation of uncertainties 
by ZEBRA such that a rescaling of these confidence in- 
terval widths by 1.25 and 1.35 for the HUDF09 and ERS 
samples, respectively, is required in order to achieve the 
68% target. Having calibrated our photometric redshift 



uncertainties we fold these via Monte Carlo simulation 
into our subsequent error calculations for the selection 
efficiencies of the color-color criteria presented in Section 

El 

In Fig. [2] we summarize graphically the key properties 
of our HUDF09 and ERS galaxy redshift catalogs. In or- 
der from the lefthand side of the page the four columns 
of this figure each contain two panels demonstrating for 
the HUDF09 and ERS (top and bottom rows, respec- 
tively): (1) the close agreement between our photomet- 
ric redshift estimates and the benchmark spectroscopic 
values where available in each sample; (11) the neces- 
sity of our (slight) re-scaling of the (ZEBRA) estimated 
uncertainties to ensure consistency in the cumulative 
distribution functions for (absolute) empirical redshift 
error divided by one sigma confidence interval width, 
i.e., I Az |/(lcr error); (ill) relative frequency histograms 
contrasting the redshift distributions of galaxies with 
spectroscopic redshifts against those for galaxies with- 
out (further subdivided according to presence/absence 
of additional GOODS-MUSIC photometry)— the latter 
distinction thereby highlighting the unique power of 
WFC3/IR imaging for exploring the z ~ 1.5 discov- 
ery space; and (iv) the overall redshift distributions for 
all galaxies in our magnitude limited samples — both of 
which peak around z ~ 1, as is typical for untargeted, 
deep imaging surveys with optical/near-IR selection (e.g. 
Rudnick et al. 2001; Ilbert et al. 2009). 

Finally, we note that our HUDF09 i^i6o < 27 mag mas- 
ter source catalog contains both z ~ 7 candidates at this 
brightness from Oesch et al. (2010) (UDFz-42566566 and 
UDFz-38807073) plus the probable supernova (UDFz- 
34537360), all of which were assigned photometric red- 
shifts z ~ 6.9 by ZEBRA. However, the remainder of 
the recently reported z ~ 7 and z ~ 8 candidates in the 
HUDF09 and ERS fields (Bouwens et al. 2010a,b; Oesch 
et al. 2010; Wilkins et al. 2010) are fainter than the ap- 
parent magnitude limits employed in this study. 

4.2. Specific Star Formation Rate and Stellar Mass 
Estimates 

Specific star formation rates (SSFRs) and stellar 
masses for all galaxies in our HUDF09 and ERS samples 
were estimated using an updated version of ZEBRA, ZE- 
BRA + (Oesch et al. 2010, in prep.). ZEBRA + computes 
various galaxy physical properties based on template 
fitting to model SEDs with the (previously-computed) 
photometric redshift, or measured spectroscopic redshift, 
held fixed. For this analysis we employed a library of 
SED templates based on the Bruzual & Chariot (2003) 
stellar population synthesis code with a Chabrier IMF 
(Chabrier 2003), drawn from a grid of exponential de- 
cay star formation rate models well-sampled over a wide 
range of ages (0.01 — 12 Gyr), decay rates (r = 0.05—9 
Gyr), and metallicities (0.05—2 Zq). The impact of 
dust reddening is accounted for during template fitting 
via a Calzetti extinction law (Calzetti 2001) with the 
E(B — V) value a free parameter (see Sargent et al. 2010 
for a demonstration of the impact of dust reddening on 
galaxies at intermediate redshift). Of course, there is a 
substantial degree of degeneracy in the choice of stellar 
population template and dust extinction model inherent 
in the fitting of broad-band SEDs, which is expected from 
previous studies to dominate the error budget (contribut- 
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Fig. 3. — The distributions of galaxies in absolute (rest-frame 
B) magnitude and stellar mass derived via the ZEBRA + code for 
galaxies in our master source catalogs to ifi60 < 27 mag in the 
HUDF09 (top row) and H\qq < 25 mag in the ERS (bottom row). 
The (conservative) SED- independent selection limits for each sam- 
ple are indicated via the green-shaded regions on each plot, while 
the boundaries of SED- dependent incompleteness (set by the differ- 
ence between the reddest and bluest spectral templates matched in 
our fits) are marked in yellow (with red hashes indicating 100% de- 
tection incompleteness). Galaxies in our samples for which reliable 
spectroscopic redshift measurements are available are highlighted 
in blue. 



ing cr\ogM ~ 0.20 dex) in samples for which broad-band 
flux measurements are available in a significant number 
of filters adequately sampling each galaxy SED (e.g. Bol- 
zonella et al. 2010; Taylor et al. 2010). Thus, for galaxies 
in our sample matched to counterparts in the GOODS- 
MUSIC survey we expect to achieve roughly this level of 
accuracy, although for galaxies with HST-only photom- 
etry our uncertainties will be necessarily higher. Details 
for ZEBRA +, including a further investigation into the 
relevant uncertainties, are given in Oesch et al. (2011, in 
prep.; see also Oesch et al. 2009). 

We note that of the 3 and 10 spectroscopically- 
confirmed passive galaxies of Daddi et al. and Cimatti 
et al., respectively, which enter our HUDF09 and ERS 
samples, all but one have ZEBRA+ estimated SSFRs 
representative of passively-evolving stellar populations 
(i.e., SSFR < 10" 10 yr _1 ). The one discrepant ob- 
ject (ERS2953) was assigned an intermediate SSFR of 
1.2 x 10" 9 yr" 1 with an extreme E(B - V) fit of 0.6 
mag. 

In Fig. [3] we present the distributions of galaxies in our 
HUDF09 and ERS master source catalogs in both abso- 
lute (rest-frame £>-band) magnitude and stellar mass as 
function of redshift. Two selection limits are indicated 
on each plot: (1) the SED- dependent completeness limit 
corresponding to the lower bound on luminosity/mass of 
the faintest /least massive systems detectable at a given 
redshift above our H\qq magnitude threshold — all mem- 
bers of our HUDF09 and ERS samples thus lie above, or 
on, this line; and (11) the SED- independent completeness 
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limit corresponding to the respective lower bound under 
the constraint of full detectability for galaxies of all spec- 
tral types. Note that the observed distributions of galax- 
ies are not tightly bound by the former until z ~ 2.5 as 
the largest K-correction from observed #160 to rest-frame 
B below this redshift occurs for the reddest spectral type 
(which is not matched to any of the faintest /least massive 
systems in these surveys; see, for comparison, Cameron 
& Driver 2007, 2009). The SED- independent limits pre- 
sented here set the thresholds (M > 1O 1O M0 for the 
HUDF09 and M > lO n M for the ERS) for selection of 
the mass-limited galaxy samples used to investigate the 
morphological mix at z ~ 2 in Section [7j 

5. OPTIMIZATION OF HST-BASED COLOR-SELECTION 
CRITERIA FOR 1.5 < z < 3.5 GALAXIES 

5.1. YHVz Color-Selection 

Through visual inspection of synthetic evolutionary 
tracks for model stellar populations in our ZEBRA+ 
template library we have identified Y105/098 — ^160 vs. 
V(306 — ^850 as the most favorable parameter space for 
selection and characterization of 1.5 < z < 3.5 galaxies 
given the set of WFC3 and ACS filters employed in the 
HUDF09/HUDF and ERS/GOODS imaging programs. 
Example tracks for instantaneous burst and constant star 
formation rate models (with metallicities of Z = 0.35 Zq 
and formation redshifts of both Zf orm = 4 and Zf orm = 6) 
are presented in Fig. [4] to illustrate the expected evolu- 
tion of active and passive systems in this plane. These 
tracks indicate that galaxies at 1.5 < z < 3.5 should 
be significantly redder in Y105/098 — #i60 , at fixed 
^606 — ^850, than galaxies at lower or higher redshifts. 
Moreover, according to the Calzetti (2001) reddening 
law, the effect of dust extinction on the observed galaxy 
colors in this diagram is simply a diagonal shift (as indi- 
cated in Fig. E| with the potential to move lower redshift 
objects into the region of high redshift color-color space 
only for passively-evolving systems at z > 1. When the 
corresponding color-color distributions of real galaxies at 
1.5 < z < 3.5 in the HUDF09 and ERS are contrasted 
against those at lower and higher redshifts (see 
again) it is clear that galaxies at these epochs are 
well separated, occupying distinct color-color sequences 
with minimal overlap. 

Guided by this qualitative analysis we define a 1.5 < 
z < 3.5 selection function, YHVz, via the combination of 
a diagonal and a vertical limit as follows: 



^105/098 - #160 > m X (V 606 ~ ^850 ) + b 
or Y] ns /oqs — #160 > c • 



105/098 



(1) 



Maximization of the selection rate (S) over this parame- 
ter space with the constraint that the z < 1 low redshift 
interloper contamination rate (C z< i) not exceed 5% re- 
turned the selection parameters m— 1.1, b — — 0.25, and 
c = 1.0 for the HUDF09, and m = 1.0, b = -0.1, and 
c = 1.5 for the ERS. Specifically, for the HUDF09: 

^105 -#160 > 1.1 x(V606-2850) -0.25 Or Fi05-#160 > 1 

(2) 
and for the ERS: 

Yb98-#160 > "^606-^850-0.1 Or Yo98-#160 > 1.5 . (3) 

These optimal YHVz selection criteria identify 96 ±1% 



and 92 ± 1% of galaxies in our target redshift interval of 
1.5 < z < 3.5 to #i6o < 27 mag in the HUDF09 and 
#i6o < 25 mag in the ERS, respectively. The corre- 
sponding strict contamination rates (C) from interlopers 
at lower and higher redshifts (z < 1.5 or z > 3.5) are only 
15 ±2% in each case, with contamination from z < 1 low 
redshift interlopers (C z< i) of only 5±1% (as constrained 
during the optimization process). The quoted uncertain- 
ties were computed according to the widths of the (me- 
dian centered) intervals containing 68% of output S, C, 
and C z< \ values recovered from a series of 1000 Monte 
Carlo simulations (with replacement) in which the col- 
ors and redshifts of galaxies were randomly shifted ac- 
cording to their estimated errors — the errors in galaxy 
colors defined according to the Poissonian errors in the 
underlying broad-band fluxes, and the redshift error s set 
to the (scaled) values from ZEBRA + (see Section 4.1) 
for galaxies with photometric redshifts only (and zero 
for those with spectroscopic redshifts). This resampling 
process also accounts naturally for the binomial errors 
(Cameron 2010) in the measured selection and contami- 
nation rates. 

5.2. Trade- Offs between Sample Completeness and 
Sample Purity 

The efficiency of YHVz selection is ultimately limited 
by intrinsic variations in the stellar population ages, star 
formation rates, and metallicities of galaxies, which pre- 
clude a perfect division by redshift on a single color- 
color diagram. Comparison of our optimal selection cri- 
teria against the model evolutionary tracks presented 
in Fig. [4] demonstrates a necessary "trade-off" between 
high completeness for both active and passive systems 
at z > 3, and exclusion of active systems at z < 1.5. 
Photometric errors in the measured colors of galaxies 
contribute only modestly to the recovered selection and 
contamination rates due to the unprecedented depth of 
the HUDF09/HUDF and ERS/GOODS imaging. Specif- 
ically, photometric errors can be excluded (above a 90% 
confidence level) as the cause of YHVz misclassifications 
for over 75% of all galaxies with conflicting spectroscopic 
redshifts in our master source catalogs, and for over 55% 
of galaxies with conflicting photometric redshifts. We ex- 
plicitly note the existence of one object in the HUDF09, 
and 3 objects in the ERS, with confident YHVz classi- 
fications 'strongly' in conflict with their measured spec- 
troscopic redshifts (i.e., offset by more than 0.5 mag to 
the wrong side of our selection limits). The one strongly 
conflicting object in the HUDF09 (UDF0045) is a com- 
pact system at z spec = 1.216 with a spectral classifica- 
tion of "broad-line AGN" (Grazian et al. 2006). Of the 
3 strongly conflicting objects in the ERS, one has an 



"AGN" spectral class (ERS2179, 



= 1.617), and 2 



are flat spectrum sources with noisy spectra for which we 
disagree with the confidence levels assigned to the red- 
shifts (ERS2585, z spec = 1.883; ERS0463, z spec = 0.117). 
The dependence of selection efficiency on apparent 
magnitude is an important consideration for the appli- 
cation of any color-selection criteria (e.g. for identify- 
ing specific high redshift galaxy types for spectroscopic 
follow up). For sub-samples of our HUDF09 and ERS 
sources at brighter magnitude limits the YHVz selection 
efficiencies remain stable but the contamination rates in- 
crease marginally. For instance, at i^i60 < 25 mag in 
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Fig. 4. — YHVz color-selection of galaxies at 1.5 < z < 3.5 in the Hubble Space Telescope ACS plus WFC3 imaging of the HUDF and 
GOODS-South fields from the HUDF09 (top row) and ERS (bottom row) programs. The observed distributions in Y105/098 — ^160 vs. 
^606 — ^850 color-color space of all objects in our HUDF09 (ifi60 < 27 mag) and ERS (ifi60 < 25 mag) master source catalogs are shown 
in the left hand panels. Galaxies with redshifts in the range 1.5 < z < 3.5 are marked with circles, while those outside this redshift interval 
are marked with triangles, and these symbols are highlighted blue and purple, respectively, if the redshift is spectroscopic, rather than 
photometric. The optimal YHVz selection criteria we derive for the high redshift population in each survey are indicated by the union of 
the gre en- a nd yellow-shaded regions — the particular significance of the yellow-shaded corner for passive galaxy selection is discussed in 
Section |5.4| Synthetic color-color tracks revealing the redshift evolution of archetypal active and passive stellar population models (with 
Z = O.35Z0, and both Zf orm = 4 and Zf orm = 6) in the relevant filters are shown for comparison in the right hand panels. Example 
reddening vectors corresponding to Calzetti (2001) extinction for galaxies observed at z = 2 are marked by black arrows. 
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Optimal YHVz 1.5 < 



TABLE 1 

< 3.5 Selection Criteria and Corresponding Selection and Contamination Rates as a Function of 
if 160 Magnitude, and Comparison Against BzK Selection^ 



Survey if i 60 Limit 



N t 



Sal 



Selection Criteria 



S (%)* C (%)* C z<1 (%) 



HUDF09 ifieo < 25 mag 309 Y 105 - H 160 > 1.1 x (Veoe - ^85o) - 0.25 or Y 105 - H i 60 > 1.0 94 ±2 17 ±4 8±2 

HUDF09 ifieo < 26 mag 540 Y 105 - H i 60 > 1.1 x (Veoe - ^85o) - 0.25 or Y 105 - H i 60 > 1.0 95 ±2 13 ±3 6±2 

HUDF09 ifi6o<27mag 986 Yi 05 - ifieo > 1.1 x (Veoe - ^sso) - 0-25 or Y 105 - H 160 > 1.0 96 ± 1 15 ± 2 5±1 

HUDF09 ifieo < 25 mag 260 (z 850 - K a ) - (B435 - z 850 ) > -0.2 or z 850 - K s > 2.5 88 ±3 34 ±3 2±2 

HUDF09 ifieo < 26 mag 395 (z 850 - K a ) - (B435 - ^sso) > -0.2 or z 850 - K a > 2.5 85 ± 2 30 ± 3 3 ± 1 

HUDF09 ifieo < 27 mag 446 (z 850 - K a ) - (B435 - z 850 ) > -0-2 or z 850 - K a > 2.5 84 ±1 30 ±2 3±1 



ERS 
ERS 
ERS 



if 160 
if 160 
if 160 



< 23 mag 

< 24 mag 

< 25 mag 



680 
1366 
2563 



Y098 — if 160 
Yogs — if 160 
Yq98 - if 160 



> V606 — ^850 ■ 

> V606 — ^850 ■ 

> V606 - ^850 ' 



0.1 or y 98 
■0.1 or y 98 
■0.1 or Yq98 



ifieo > 1-5 
ifieo > 1-5 
ifieo > 1-5 



88 ±4 
92 ±2 
92 ±1 



21 ±4 
15 ±2 
15 ±2 



10 ±3 
5±1 

5±1 



ERS ifieo < 23 mag 520 (z 850 - K B ) - (B435 ~ z 850 ) > -0.2 or z 850 - K s > 2.5 84 ±4 38 ±3 6±3 

ERS ifieo < 24 mag 1009 (z 850 - K B ) - (B435 - z 850 ) > -0.2 or z 850 - K s > 2.5 90 ± 2 31 ± 2 3 ± 1 

ERS ifieo < 25 mag 1768 (z 850 - K a ) - (B435 - z 850 ) > -0.2 or z 850 - K s > 2.5 86 ± 1 30 ± 2 4±1 



tOnly galaxies matched to counterparts in the GOODS-MUSIC catalog with measured K a magnitudes are used in the analysis of the BzK 
selection efficiency. 

*For YHVz selection a target redshift interval of 1.5 < z < 3.5 is used for computing the selection efficiencies and strict contamination 
rates, while for BzK selection we suppose a target redshift interval of 1.4 < z < 2.5. 



the HUDF09 and H 160 < 23 mag in the ERS we re- 
cover contamination rates of C = 17 ± 4% and 21 ±4%, 
and C 2< i = 8 ± 2% and 10 ± 3%, respectively. This 
is an expected consequence of the higher ratio of low 
redshift sources to high redshift sources in the brighter 
sub-samples, given the intrinsic limitations of redshift de- 
termination in a single color-color diagram. The depen- 
dence of YHVz selection efficiency on apparent magni- 
tude in our HUDF09 and ERS source catalogs is further 
quantified in Table [I] 

5.3. A Comparison of Performance 

We have compared the efficiency of YHVz 1.5 < z < 
3.5 selection against the BzK 1.4 < z < 2.5 selection 
of Daddi et al. (2004). To this purpose, we have com- 
puted the corresponding BzK selection and contamina- 
tion rates for all galaxies in our master source catalogs 
well-matched to GOODS-MUSIC objects with measured 
K s magnitudes. For BzK we obtain S = 84 ± 1% and 
86±1%, C = 30±2% and 30±2%, and C z<1 = 3±1% and 
4 ± 1% to #i 6 o < 27 mag in the HUDF09 and H 160 < 25 
mag in the ERS, respectively (as summarized in Table 
HI). Comparatively, YHVz selection offers a more effi- 
cient mechanism for high redshift galaxy selection with 
a significantly more precise sampling of the target red- 
shift interval than BzK (i.e., lower strict contamination 
rates) . 

A further distinction between YHVz and BzK is that 
the former fully exploits the very high angular resolution 
and uniformly high sensitivity of both the WFC3 and 
ACS, whereas the latter requires K s imaging, not avail- 
able, for large galaxy samples, at a similarly high resolu- 
tion. As a result, only YHVz selection allows the confi- 
dent photometric identification of high redshift galaxies 
with neighbors at angular separations below the confu- 
sion limit of ground-based K s -band imaging. This en- 
ables more detailed studies of the structural assembly 
and star formation properties of galaxies at these crucial 
epochs. 

5.4. Passive YHVz 

According to the stellar population evolution models 
presented in Fig. [5j galaxies at 1.5 < z < 3.5 with ac- 
tively star-forming and passively-evolving stellar popu- 
lations occupy distinct regions of the Y105/098 — ^160 vs. 
Veoe — ^850 color-color plane. In particular, high redshift 
passively-evolving systems are expected to exhibit much 
redder Vqoq — zs5o and Y105/098 — ^160 colors than star- 
forming systems at the same epochs. This is consistent 
with the observed colors of spectroscopically-confirmed 
passive galaxies at 1.5 < z < 3.5 from the Daddi et al. 
and Cimatti et al. samples, as shown in Fig. [5] Inter- 
estingly, all but one of these confirmed passive galaxies 
can be identified simply via an extension of our optimal 
YHVz selection criteria to define a passive YHVz selec- 
tion analagous to the pBzK object class. Specifically, for 
theMJDF09: 



^105-^160 < l.lx(V606-285o)-0.25 and Fi05-#160 > 1 

(4) 
and for the ERS: 



Yr 



098 



#160 < Vcso6 ~ ^850 - 0.1 and F 98 - #ieo > 1-5 . 

(5) 

In addition to the well-studied pBzK objects, our pas- 
sive YHVz selection criteria identify an additional 4 and 
10 galaxies with photometric, and/or spectroscopic, red- 
shifts of 1.5 < z < 3.5 in the HUDF09 and ERS, re- 
spectively. (The one passive YHVz galaxy with a spec- 
troscopic redshift not from either the considered pBzK 
catalogs, ERS2182, exhibits only a tentative Oil emis- 
sion line detection at the upper wavelength limit of its 
VLT/FORS2 optical spectrum from the GOODS pro- 
gram; Vanzella et al. 2005, 2006, 2008). Photometric er- 
rors in the colors of these additional passive candidates 
enable us to robustly assign 8 (4 from the HUDF09, and 
4 from the ERS) to the passive YHVz class above a 90% 
confidence level. 

According to the Calzetti (2001) model, dust reddening 
should move galaxies at z ~ 2 near-parallel to our passive 
YHVz boundary, implying minimal contamination of our 
passive sample from dust reddened, star-forming objects 
at 1.5 < z < 3.5. Indeed, the best-fit model templates to 
the observed SEDs of all confident passive YHVz galax- 
ies at 1.5 < z < 3.5 in our catalogs were representative 
of very low specific star formation rates (SSFR <10 -10 
yr -1 ), with only one galaxy requiring a high dust ex- 
tinction value (E(B — V) > 0.3 mag) for an optimum 
fit. Only 3 of the 10 additional passive YHVz candi- 
dates in the ERS not in either pBzK sample have ZE- 
BRA + SEDs favouring models with intermediate SSFRs 
(10" 10 < SSFR <10" 9 yr -1 ) and high dust reddening val- 
ues; all 3 enter our selection below the 90% confidence 
level. Similarly, 4 of the 17 total (pBzK plus new) passive 
YHVz systems in the ERS were identified as strong MIPS 
24 /im sources (F 24/im > 20 /xJy) in the GOODS-MUSIC 
catalog (see Fig.|5|. The detection of strong 24 /im emis- 
sion from galaxies at these epochs provides an alternative 
indicator of the presence of dust-obscured star- format ion 
(cf. Yan et al. 2004), modulo any contribution from AGN 
activity; again, of the four 24 /im sources contaminating 
our passive YHVz catalogue, three make our selection 
only at low confidence levels. 

In Fig. [6] we present WFC3 i^i6o postage stamp im- 
ages of all passive YHVz galaxies at 1.5 < z < 3.5 from 
our HUDF09 and ERS master source catalogs. Visual 
inspection confirms that these are typically compact sys- 
tems with highly centrally-concentrated light profiles and 
smooth morphologies, in agreement with previous studies 
of this galaxy class (cf. Van Dokkum et al. 2004, 2010; 
Daddi et al. 2005; Trujillo et al. 2006; Cimatti et al. 
2008; Franx et al. 2008; Toft et al. 2009; Damjanov et al. 
2009; Smozoru et al. 2010). However, at least 6 of these 
passive YHVz galaxies exhibit morphologies suggestive 
of the presence of either a stellar disk with or without 
clumpy sub-structure (UDF0681, ERS0177, ERS1000, 
UDF0093, ERS3054). The presence of disk-like mor- 
phologies within a subset of passive z ~ 2 galaxies has 
previously been reported by Stockton et al. (2008), Mc- 
Grath et al. (2008), and Cassata et al. (2010). 

We stress that the passive YHVz selection criteria (sim- 
ilarly to, for example, the pBzK criterion), may not re- 
turn the complete sample of passively evolving galax- 
ies at the targeted redshifts. At face value, there are 
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Fig. 5. — Passive YHVz color-selection of passively-evolving galaxies at 1.5 < z < 3.5 in the Hubble Space Telescope ACS plus WFC3 
imaging of the HUDF and GOODS-South fields from the HUDF09 (left column) and ERS (right column) programs. The Y 105 / 098 — ^160 vs. 
V606 — ^850 color-color distributions of galaxies at these redshifts with SED-fit classified actively star-forming stellar populations are marked 
with blue circles and passively-evolving stellar populations with purple triangles. Galaxies with SED-fits requiring high dust reddening 
values (E(B — V) > 0.3 mag) are marked with small white dots, and strong MIPS 24 /im sources (F^A^m > 20 /iJy) are marked with 
black circles. Yellow shading is used to highlight the sub-region of color-color space defined by our passive YHVz selection criteria. The 
photometric errors in the measured colors of 1.5 < z < 3.5 galaxies within this region are illustrated by the corresponding error bars, and 
the positions of spectroscopically-confirmed passive galaxies at z > 1.5 from the pBzK catalogs are highlighted with red circles. 



galaxies with very low SSFRs computed with ZEBRA + 
SED fits, in both the HUDF09 and ERS samples at 
1.5 < z < 3.5, which reside outside our passive YHVz 
limits, occupying a region of intermediate color-color 
space shared with dusty, star- forming sources. Given 
that this 'green valley' region is consistent with the im- 
pact of moderate-to-high dust reddening on star-forming 
systems at 1.5 < z < 3.5, however, we suspect that 
the inherent degeneracy between the choice of reddening 
value and stellar population template during SED-fitting 
is responsible for the under-estimation of SSFRs in many 
of these additional passive candidates. 

5.5. Other HST-based Criteria for Selecting z > 1.5 
Galaxies: 
The JHVI Criterion 

Motivated by the desire to capitalize on other HST 
high-z surveys employing similar but not identical pass- 
bands, we have explored the performance of other color 
selection criteria in isolating galaxies in the z > 1.5 red- 
shift regime. In particular, the CANDELS team is cur- 
rently in the process of surveying a total area of ~800 
arcmin 2 on the sky with WFC3/IR and ACS with a two- 
tier observing strategy featuring a 'wide, but shallow' 
component (CANDELS/Wide) and a 'deep, but narrow' 
component (CANDELS/Deep; as described in Grogin et 
al. 2011 and Koekemoer et al. 2011). CANDELS/Deep 
features observations of regions in the two GOODS fields 
in each of the Y105, ^125, and Hiqq filters, supplement- 
ing the existing GOODS £^35, V<306j ^775 1 an d ^850 ACS 
images. Hence, our YHVz selection can be readily ap- 
plied for the study of 1.5 < z < 3.5 galaxies in this 
deep dataset. CANDELS/Wide, however, will feature 



shallower imaging of regions from the GOODS-N & S, 
COSMOS, EGS, and UKIDSS/UDS fields in only the 
^606 ? ^814, ^125, and #160 filters, precluding application 
of our YHVz selection criteria (except in the GOODS-N 
region for which Y105 data is available from the ERS and 
CANDELS/Deep programs). 

To estimate the efficiency of other potential color- 
color selection strategies for CANDELS/Wide, we first 
estimate isi4 magnitudes for all objects in our ERS 
Hiqq < 25 mag source catalog by interpolating the best- 
fit SED template for each galaxy against the 2gi4 trans- 
mission function. The wavelength coverage of the igi4 fil- 
ter is significantly longer and redder than that of the ^775 
filter employed in the GOODS ACS program, providing 
deeper imaging in this part of the spectrum for a given 
exposure time. Using these interpolated Jsi4 magnitudes 
we contrast the distributions of galaxies at z > 1.5 and 
z < 1.5 in our ERS #160 < 25 mag sample in J125 — H\qo 
vs. ^606 — ^814 color-color space in Fig. [7] 

By re peat ing the optimization procedure employed in 
Section |5.1| we recover the following JHVI selection cri- 
teria: 

J125-H16O > 1.5x(V606-^8i4)-0.1 or J125-H16O > 0.5 . 

(6) 
It is clear from visual inspection that the distinction 
between galaxies above and below z ~ 1.5 in this JHVI 
parameter space is less efficient than that demonstrated 
earlier for our YHVz selection. Nonetheless, these crite- 
ria do identify 77 ± 2% of ERS #160 < 25 mag galaxies 
at z > 1.5 with a strict (z < 1.5) contamination rate 
of 26 ± 1% and aKl interloper rate of only 5 ± 1% 
(as constrained in the optimization). Thus, despite the 
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Fig. 6. — Hieo postage stamp images of passive YHVz galaxies at 1.5 < z < 3.5 in the HUDF09 (top row) and ERS (bottom three rows) 
WFC3 imaging of the HUDF and GOODS-South fields. In each case the elliptical Kron aperture identified for that galaxy is overlaid on 
the image, with the ellipse color revealing whether its best-fit SED template was indicative of a passively-evolving (purple) or star-forming 
(blue) stellar population. Solid ellipses are used to denote galaxies classified within the passive YHVz selection area at >90% confidence 
given the photometric errors on their colors, while broken ellipses are used for lower confidence classifications. The corresponding H\qq 
apparent magnitude and best available redshift for each source are noted in the top and bottom right-side corners of each postage stamp 
image, respectively, with spectroscopic redshifts highlighted in yellow. Additional notes are added for sources spectroscopically-confirmed 
as passively-evolving systems in the pBzK Cimatti et al. or Daddi et al. samples, or alternatively, for sources identified as strong MIPS 24 
/im detections (F 24/ xm > 20 /iJy) in the GOODS-MUSIC catalog. 



lower performance relative to our main YHVz selection 
criteria, this alternative JHVI selection still provides a 
useful tool for the identification of z > 1.5 galaxies. 

5.6. Lowering the Low-z Limit: 
VIH Galaxies Above and Below z = 1 

As a final remark, we note that by modifying the goals 
of the selection to z > 1 galaxies one can identify more 
efficient color-color criteria using only the V606 5 ^814, and 
Hiqq filters, as shown in Fig. [7[ The optimum VIH cri- 
teria we recover for the ERS source catalog are: 



hu ~ #160 > 1-75 x (V 606 - Jsu) - 0.25 
or J 814 -#i6o > 2.1 . 



(7) 



These criteria detect 85 ± 1% of ERS #160 < 25 mag 
galaxies at z > 1 (with only 5 =b 1% contamination from 
z < 1 interlopers). 

6. SOURCE CATALOGS 

The complete source catalogs employed in this study, 
containing all robust detections to #160 < 27 mag in the 
HUDF09 and H 160 < 25 mag in the ERS, are publicly- 
available in electronic form from the following URL: 



http : //www . astro . phys . ethz . ch/cameronetal201 1/ 
. These catalogs have been carefully inspected both to 
remove spurious detections and to ensure selection of the 
optimum photometric aperture for flux measurement of 
each object. A series of flags is used to indicate which 
objects appear within the area covered by the GOODS- 
MUSIC catalog, which are classified as well-matched to 
a counterpart in the GOODS-MUSIC catalog, which 
appear in a region of high signal-to-noise imaging in all 
7 ACS plus WFC3 filters, which have a GOODS-MUSIC 
catalog spectroscopic redshift significantly in conflict 
with their observed colors, which are spectroscopically- 
confirmed stars or unresolved stellar candidates, and 
which are in our passive YHVz selection. The column 
numbers and definitions for these catalogs are listed in 
Table |2] below. 

7. COMPARISON OF REST-FRAME OPTICAL VERSUS 
REST-FRAME UV MORPHOLOGIES OF 1.5 < z < 3.5 
GALAXIES 

As a first application of our i^i6o-selected HUDF09 
and ERS source catalogues presented in Section [6] we es- 
tablish the most basic, and thus fundamental, facts con- 
cerning the morphologies of 1.5 < z < 3.5 galaxies in the 
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Fig. 7. — Potential strategies for color-selection of high redshift galaxies in combinations of WFC3 and ACS filters from the CAN- 
DELS/ Wide program (excluding Yios which is available only for the GOODS-N field). Specifically, we demonstrate the separation of 
galaxies at z > 1.5 (circles) and z < 1.5 (triangles) in J125 — iii60 vs. V606 — isi4 color-color space (left hand panel), and at z > 1 (circles) 
and z < 1 (triangles) in 7si4 — iii60 vs. V606 — isi4 color-color space (right hand panel), using our ERS if 160 < 25 mag sample. In each 
case, colored symboles are used if the redshift is spectroscopic, rather than photometric. The optimal JHVI and VIH selection criteria 
we derive for the hi gh re dshift population in each filter combination are indicated by the union of the green- and yellow-shaded areas. As 
described in Sectio n|5.5| the isi4 filter magnitudes of these ERS sources were estimated by interpolating the best-fit SED template for each 
galaxy against the isi4 transmission function. 



rest-frame optical window compared with the morpholo- 
gies of the same galaxies in rest-frame UV light. 

In recent years, ground-based or HST/NICMOS NIR 
imaging sampling rest-frame optical wavelengths at 
1-5 ;$ z < 3.5 has uncovered a galaxy population 
at these epochs rich in diversity of structure: with 
smooth disk/bulge-disk galaxies, and rather normal ellip- 
ticals coexisting with massive interacting/disturbed sys- 
tems, as well as ultra-compact, elliptical-like galaxies, 
large, turbulent, clumpy disks, and compact, dispersion- 
dominated disks (Corbin et al. 2000; Daddi et al. 2005; 
Papovich et al. 2005; Franx et al. 2008; Stockton et al. 
2008; Genzel et al. 2008; Conselice et al. 2008; Buitrago 
et al. 2008; Forster-Schreiber et al. 2009; Elmegreen et 
al. 2009). Interestingly, many galaxies at z > 1.5 exhibit 
genuinely disturbed/interacting morphologies, i.e., not 
only the young stellar populations visible in UV light, 
but, rather, most of the stellar mass in these galax- 
ies, is arranged in an irregular manner. With the sig- 
nificantly higher angular resolution and uniformly high 
imaging sensitivity of WFC3/IR, the new HUDF09 and 
ERS #160 images offer the most detailed view to-date of 
the high-redshift galaxy population at rest-frame optical 
wavelengths. Through visual classification of the galax- 
ies in our 1.5 < z < 3.5 color-selected HUDF09 and ERS 
samples, we thus investigate here, qualitatively, the di- 
versity of galaxy structure as revealed by the WFC3/IR 
#i6o (rest-frame optical) images relative to the ACS 2:350 
(rest-frame UV) images. 

The first indication of the paramount importance of 
the rest-frame WFC3 images in understanding the 1.5 < 
z < 3.5 galaxy population comes from the comparison 
with the GOODS-MUSIC catalog discussed in Section 



3.1 In most cases where no counterpart to one of our 
#160 sources was found in the GOODS-MUSIC cata- 
log, the centroid offset that prevented automatic catalog 
matching was due to morphological differences between 
the rest-frame UV zs^o image (used as a basis for most 
GOODS-MUSIC detections) and the new rest-frame op- 
tical #160 image provided by the WFC3. In a number 
of instances, the greater sensitivity of the WFC3 #i60 
image to the smoother, underlying stellar mass distri- 
bution provided a confident aperture definition and cen- 
troid conflicting significantly with that adopted in the 
GOODS-MUSIC analysis, and no match was possible. It 
is evident that the rest-frame UV images on which the 
GOODS-MUSIC catalog was based often identify, as the 
galaxy center, a knot of star formation substantially dis- 
placed from the true galaxy center. This is, on the other 
hand, well identified in the rest-frame optical WFC3 im- 
ages. 

Also, the combination of passive- SED color-selection 
(passive YH Vz) and rest-frame optical morphologies (see 
Section |5.4| reveals that passive galaxies at such early 
epochs are indeed mostly spheroids, as shown in previous 
studies (e.g. Van Dokkum et al. 2004, 2010; Daddi et al. 
2005; Trujillo et al. 2006; Cimatti et al. 2008; Franx et al. 
2008), but, in addition, a substantial fraction (~15-20%) 
show disk morphologies. We speculate that this popula- 
tion of 'quenched' disks at such early epochs might be a 
smoking gun for disk galaxies which have been cut out 
of their feeding streams of cold gas, as seen in numeri- 
cal simulations and expected on theoretical grounds (e.g. 
Keres et al. 2005; Dekel et al. 2008). Mergers of such 
disks would be 'dry' while nonetheless involving galaxies 
with a dynamically-cold disk component. Foreseen large 
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Fig. 8. — Rest-frame UV (ACS Z850) and optical (WFC3 if 160) postage stamp images of example 1.5 < z < 3.5 galaxies revealed as disk, 
or bulge/disk, systems by WFC3 (CLASS 1). In each case the elliptical Kron aperture identified for that galaxy is overlaid on the image. 
The corresponding Z850 or ifi60 apparent magnitude and best available redshift for each source are noted in the top and bottom right-side 
corners of each postage stamp image, respectively, with spectroscopic redshifts highlighted in yellow. 




Fig. 9. — Same as for Figure[5] but this time for CLASS 2, 1.5 < z < 3.5 example galaxies that show a disk-less, spheroid-only morphology 
both in the rest-frame optical WFC3 images and in the rest-frame UV ACS images. 



i7i6o-band surveys with the HST will be instrumental in 
understanding what fraction of dry mergers at z ~ 2-3 
involve disk galaxies — an important constraint for galaxy 
evolution models. 

Furthermore, by comparison of the ACS and WFC3 
morphologies of all galaxies in our master source catalogs 
within the 1.5 < z < 3.5 redshift window we identify four 
main galaxy sub- types. The first category (CLASS 1) are 
galaxies which, in the i^i6o~band WFC3 images, reveal 
a regular disk morphology (with or without a bulge) — in 
contrast with their irregular rest-frame UV morpholo- 
gies. These CLASS 1 galaxies are the only ones which 
show a morphological difference between their rest-frame 
UV and rest-frame optical images. The remaining three 
classes do separate galaxies with very different morpholo- 
gies, but, strikingly, these morphologies are unchanged in 
rest-frame UV and rest-frame optical light. These classes 
are, respectively: galaxies which, both in the rest-frame 
UV and rest-frame optical NIR WFC3 images maintain 
a regular, spheroid-like morphology (CLASS 2); galaxies 
which retain, in the NIR WFC3 images, a highly irreg- 
ular/disturbed morphology, similar to that observed in 
the ACS rest-frame UV images (CLASS 3); and, finally, 
the population of small-sized, smooth and rather sym- 



metric 'blobs' of light in both the #160 and zs5o images 
for which it is unclear whether these systems are proto- 
disks or spheroids (CLASS 4). Examples of each class are 
given in Figures [§J J9j [lOj and [TT] 

In CLASS 1 objects, the rest-frame optical morphol- 
ogy is substantially smoother, i.e., less clumpy, and 
also more centrally concentrated, than in the rest-frame 
UV. Indeed, in such systems, the WFC3 i^i6o im- 
ages expose a well-defined center — in contrast with the 
clumpy/asymmetric appearance of these galaxies in the 
rest-frame UV light. Two such galaxies in our sam- 
ple (UDF0093, ERS2047) exhibit clear bulge-dominated 
morphologies, and a number of others (e.g. UDF0590, 
ERS3032) indicate the presence of small, central, disk 
light excesses, most likely small (er) bulges, consistent 
with counterparts identified by Papovich et al. (2005) 
and Stockton et al. (2008). 

Class 2 objects show the smooth and compact ap- 
pearance indicated by previous studies to be typical for 
many spheroids at these epochs (e.g. Daddi et al. 2004, 
2005; Cimatti et al. 2008; Franx et al. 2008; Damjanov et 
al. 2009). Already Szomoru et al. (2010) have used the 
HUDF09 WFC3/IR images to establish, on firm grounds, 
the compactness of one 'quenched' spheroid at z = 1.98 
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Fig. 10. — Same as for Figure|5] but this time for CLASS 3, 1.5 < z < 3.5 example galaxies that keep a similar, highly irregular/disturbed 
morphology in rest-frame opticaTas in the rest-frame UV. 




Fig. 11. — Same as for Figure [5] but this time for CLASS 4, 1.5 < z < 3.5 example galaxies that have a smooth, 'blob' morphology both 
in the WFC3 if 160 an d ACS Z850 images. 



(UDF0884), extracted from the Daddi et al. sample. 
While many galaxies in this morphological class are iden- 
tified as passive YH Vz s ystems in our color-color selec- 
tion criteria (Section |5.4| ), we stress that a non-negligible 
fraction of them have instead star- forming SEDs, indicat- 
ing that such smooth, compact morphologies, while more 
typical of the passively-evolving galaxy population at 
these redshifts, can also, to a smaller extent, be achieved 
by galaxies that are still actively forming stars at those 
epochs (e.g. ERS0556, ERS1202). 

Class 3 objects exhibit a wide variety of morphologies, 
from systems with a compact nucleus plus a highly asym- 
metric disk-like structure (e.g. UDF0124) to systems 
with multiple bright nuclei (e.g. UDF0435, ERS2010, 
ERS1203). Interestingly, the disk-like structures in these 
systems are often warped (e.g. UDF0124, UDF0870), a 
phenomenon that is observed to occur in high-resolution, 
hydrodynamical numerical simulations in CDM cosmol- 
ogy, as disk galaxies stream down the hot gas filaments of 
the cosmic web (Steger et al 2010, in prep.). In a number 
of cases (e.g. UDF0435, UDF0870) there exist clear dif- 
ferences between the ^50 and #160 brightnesses of spe- 
cific nuclei within these galaxies, reflecting a complex, 
spatially-dependent — i.e., 'clumpy' — star formation ac- 
tivity, and thus history, in these systems. Since these 
systems will most likely evolve into relatively massive 



galaxies by redshift zero, a substantial stellar population 
mixing, through mergers or internal dynamical instabili- 
ties, must occur over cosmic time to achieve the massive 
galaxy population with rather smooth stellar population 
gradients that we observe today. 

One CLASS 3 irregular galaxy in our sample, a lop-sided 
disk (UDF0870) at z = 1.509, might exhibit a large-scale 
bar. While a decline in the bar fraction with increasing 
redshift has been demonstrated by several studies out to 
redshifts of z ~ 1 (e.g. Abraham et al. 1999; Sheth et 
al. 2008; Cameron et al. 2010), it is only with the newly 
available rest-frame optical, deep and high-resolution im- 
ages of high-z galaxies that a firm statement on the bar 
fraction at z > 1 will become feasible. The discovery of 
one such system in the HUDF09 images indicates that 
large-scale, coherent dynamical instabilities such as bars 
can and do exist at such early epochs (see also the z > 2.5 
'bar-like structures' classified in the rest-frame UV by 
Ravindranath et al. 2006), although they may be not as 
frequent as at redshift z ~ 1, given that disks at such 
early times appear to be, on average, substantially hot- 
ter and more turbulent than at later times (e.g. Genzel 
et al. 2006, 2008, 2010). 

Class 4 objects are relatively small-sized systems ex- 
hibiting smooth, symmetric morphologies, but with less 
centrally-concentrated light distributions than the CLASS 
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TABLE 2 
Column Numbers and Definitions for HUDF09 and ERS Source Catalogs 



# Definition 



1 ID (in format UDF#### / ERS####) 

2 Right ascension (deg) (J2000) 

3 Declination (deg) (J2000) 

4 if 160 Kron magnitude (AB)t 

5 Hiqo RMS error in Kron magnitude"^ 

6 if 1 60 Kron radius (arcsec)t 

7 if i 6 o Ellipticity (l-b/a) + 

8 if 160 Position angle (deg) (for images aligned such that west=0 deg and north=90 deg)t 

9 if 160 Stellarity parameter^ 

10 B435 Kron magnitude (AB)t 

11 S435 RMS error in Kron magnitude (99 if null detection) ^ 

12 V606 Kron magnitude (AB)t 

13 Vqoq RMS error in Kron magnitude (99 if null detection) ^ 

14 f775 Kron magnitude (AB)t 

15 f775 RMS error in Kron magnitude (99 if null detection) ^ 

16 ^850 Kron magnitude (AB)t 

17 ^850 RMS error in Kron magnitude (99 if null detection) t 

18 I105/098 Kron magnitude (AB) (Y 105 in HUDF09, Yo98 in ERS) 1 " 

19 ^105/098 RMS error in Kron magnitude (99 if null detection) t 

20 J125 Kron magnitude (AB)t 

21 J125 RMS error in Kron magnitude (99 if null detection)''' 

22 GOODS-MUSIC match flag (0=unmatched, l=matched) 

23 GOODS-MUSIC ID (0 if unmatched) 

24 GOODS-MUSIC area flag (0=not in area, l=in area) 

25 Object appears in a high signal-to-noise region in all 7 bands flag (0=not in high S/N region, l=in high S/N region)* 

26 Flag if GOODS-MUSIC catalog spectroscopic redshift is strongly in conflict with the observed color (0=okay, l=in conflict) 

27 Star-galaxy separation flag (0=galaxy, Inspect roscopically-confirmed star, 2=candidate star) 

28 ZEBRA photometric redshift 

29 x 2 °f photometric redshift template fit 

30 Number of filters used in photometric redshift template fit 

31 Passive YHVz selection flag (0=not in selection, l=in selection) 

t Parameters in columns 4-21 derived from SExtractor output; see Bertin &; Arnouts (1996) for detailed definitions. 

* Regions of "high signal-to-noise" imaging defined as the area on each frame with local median RMS less than twice the global median 

RMS. 
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Fig. 12. — Fractional contribution of high- z morphological sub- 
types as a function of stellar mass at 1.5 < z < 2.15 (top row) and 
2.15 < z < 3.5 (bottom row) from visual inspection of the ACS 
z 8 50 and WFC3 H 160 images of 30 galaxies at M > 10 10 M o in 
the HUDF09 and 30 galaxies at M > 10 n M o in the ERS. Green, 
yellow, blue, and purple histograms indicate, respectively, disks 
unveiled by WFC3 (CLASS 1), spheroids (CLASS 2), highly irregu- 
lar/disturbed UV and optical morphology (CLASS 3), and small- 
sized 'blob'-like systems (CLASS 4). 

2 objects. As such, these CLASS 4 objects could either 
be diffuse counterparts of the CLASS 2 spheroids, or else 
small proto-disks in the earliest stages of growth. Future 
studies of the quantitative structural parameters and, 
ideally, kinematical properties of this object class are re- 
quired to establish the place of these systems within the 
cosmic history of galaxy structural assembly. 

In Fig. [12] we present the morphological mix of these 
galaxy types in each of the ~1.2 Gyr intervals, 1.5 < 
z < 2.15 and 2.15 < z < 3.5, based upon counts in 
mass-limited samples of 30 galaxies at M > 1O 1O M 
in the HUDF09 and 30 galaxie^p] at M > 10 n M o in 
the ERS (see Fig. [3] for the relationship between our 
mass completeness limits and the distribution of galax- 
ies in our HUDF09 and ERS source catalogs as a func- 
tion of redshift). Interestingly, as we confirm via arti- 
ficial noise degradation analysis of the HUDF09 images 
in the Appendix to this paper, despite a ~14x shorter 
exposure time the depth of the ERS imaging is sufficient 
to perform morphological classification of the (relatively 
crude) nature described here at comparable accuracy to 
that in the HUDF09. The combined (i.e., >1O 1O M ) 
fractions of objects in the four different classes are, re- 
spectively, 29±?%, 45±|%, 13±|% and 13±|% at 1.5 < 
z < 2.15, and 0±g%, 18±§%, 41±|% and 41±|% at 
2.15 < z < 3.5 (with la binomial uncertainties estimated 
using a Bayesian approach; cf. Cameron 2010). Note that 
the evolution of these fractions likely implies substantial 
morphological transformations across the redshift z ~ 2 
epoch of peak star formation activity. We highlight, in 

10 Note that it is entirely by coincidence, not design, that there 
exist exactly 30 galaxies in each of these mass-limited samples! 



particular, the emergence of a substantial population of 
regular disk morphologies, and a substantial increase in 
the fraction of spheroidal morphologies at epochs z < 2. 

Importantly, the four morphological classes show dif- 
ferent redshift and stellar mass distributions. At 2.15 < 
z < 3.5 CLASS 3 objects (with highly irregular mor- 
phologies also in rest-frame optical light) and CLASS 4 
objects (with small-sized, 'blob'-like morphologies) are 
overwhelmingly dominant at both intermediate and high 
mass scales, while CLASS 2 spheroids are largely re- 
stricted to above 10 11 M . At 1.5 < z < 2.15 CLASS 3 
and CLASS 4 objects are now largely restricted to inter- 
mediate masses, while the CLASS 2 spheroid population 
now dominates at high masses. At 1.5 < z < 2.15 the 
Class 1 systems emerge and contribute ~25% of the 
morphological mix at both mass scales. At and around 
the epoch of peak cosmic star formation rate in the uni- 
verse, therefore, the statistical census of the brightest 
galaxies indicates that the most massive of them are as- 
sembled in regular, Hubble- type elliptical and disk (plus 
bulge) rest-frame optical morphologies, while less mas- 
sive systems often reveal, in their rest-frame optical im- 
ages, clear signs of dynamical youth. 

Perhaps the most striking morphological difference be- 
tween these two redshift intervals is the transition across 
the z ~ 2 epoch for galaxies at the highest stellar masses, 
>1O 11 M , from a diversity of structural types to a dom- 
inant population of spheroids. Although SED-fit based 
SSFRs suffer from large uncertainties, and our samples 
of galaxies with masses above lO n M at these red- 
shifts are only small (16 at 2.25 < z < 3.5 and 15 at 
1.5 < z < 2.25), we note that the assembly of the massive 
spheroids appears to be accompanied by a rapid 'quench- 
ing' of star formation. Adopting a redshift-dependent di- 
vision between star-forming and passively-evolving sys- 
tems based on the SSFR necessary for a galaxy to dou- 
ble its stellar mass between its observed redshift and the 
present day, and employing o ur SE D-fit SSFRs computed 
with ZEBRA + (see Section 4.2), we find the fraction 
of star- forming galaxies at >lO ii M falls from 63±^3% 
at 2.25 < z < 3.5 to b±f at 1.5 < z < 2.25. Us- 
ing our own color-based criteria for identification of a 
"passively-evolving" SED, i.e., passive YHVz, the frac- 
tion of quenched galaxies above 10 10 M increases from 
6±f% at 2.25 < z < 3.5 to 60±^% at 1.5 < z < 2.25. 
This is consistent with the observed decline in the star- 
forming fraction estimated by Feulner et al. (2005) for 
a selection of galaxies in the FORS-Deep and GOODS- 
South fields. If confirmed by larger samples, these re- 
sults would indicate that z ~ 2 is the epoch when both 
morphological transformations first produce a significant 
population of massive spheroids and physical processes 
quench their star formation activity and move them to 
the red sequence, as expected in recent models (Peng et 
al. 2010). 

8. SUMMARY 

The wealth of WFC3 imaging of the high redshift 
universe achieved or planned, coupled with archival 
ACS optical imaging, is self-sufficient, without additional 
ground-based photometry, to robustly identify galaxies 
in the 1.5 < z < 3.5 redshift window. In this paper 
we have presented, to this aim, new, robust and effi- 
cient color-selection criteria established primarily on the 



17 



^105/098 - ^160 vs. Veoe - ^850 color-color diagram (and 
similar ones, slightly less efficient but still valuable, which 
use combinations of the HST J125, #160? ^606 5 an d hi4 
filters, or simply the V606, ^814, and #160 filters). We 
have also shown that an additional color selection on the 
^105/098 - #160 vs. V606 - ^850 color-color diagram fur- 
ther enables the disentanglement of passively-evolving, 
'quenched' galaxies from actively star-forming galaxies 
at these redshifts. 

As a first use of the new galaxy catalog that we have 
produced, starting from the i^i6o WFC3 images, we have 
performed a qualitative comparison of 1.5 < z < 3.5 
galaxy morphologies between rest-frame UV and rest- 
frame optical light. Regular disk galaxies — with or 
without a bulge — are present at these epochs — the disk 
components being unseen or irregular in rest-frame UV 
images, and being clearly revealed in the rest-frame 
optical images. Together with the population of smooth 
spheroids, these regular disk galaxies have substantial 
stellar masses, >1O 1O M (or, for the spheroids, even 
above this threshold). At masses >10 n M o , the WFC3 



NIR images reveal a transition of morphologies from a 
variety of structures at z > 2 to a marked dominance 
of elliptical-like stellar distributions at z < 2. Across 
the z ~ 2 boundary the fraction of star-forming galaxies 
at >1O 11 M decreases from ~60% to ~5%. Consistent 
with analysis of the z < 1 population (Peng et al. 2010), 
and observed here in rest-frame optical morphologies 
and colors, z ~ 2 appears therefore to be the epoch at 
which assembly and quenching of the massive spheroid 
population begins in earnest. The WFC3 NIR images 
also reveal, however, that, less massive galaxies at the 
same epochs often keep a highly irregular morphology, 
also in the rest-frame optical light, as they had in the 
rest-frame UV light. This indicates that their previously 
reported irregular morphologies in rest-frame UV are 
not driven by detection of star-forming populations on 
top of settled, underlying older stellar populations, but 
rather an indication of a genuine juvenile dynamical 
state for much of the 10 10 < M < 10 n M© galaxy 
population at these redshifts. 
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APPENDIX 

In this brief appendix we confirm through statistical image degradation that the limited depth of the ERS imaging 
with WFC3/IR is indeed sufficient to allow for confident morphological classification of z ~ 2 galaxies in the rest-frame 
optical, at least to the relatively crude level of detail employed in the analysis of Section [7[ Adopting the ultra-deep 
HUDF09 #160 images of our 1.5 < z < 3.5 master sample as the benchmark for morphological classification at these 
redshifts we thus aim to demonstrate a consistency of visual appearance following artificial degradation to the noise 
level of the ERS. In particular, the net exposure time of the master UDF #160 image is a factor of ~14 times greater 
than that of the ERS #i6o i ma g e 5 indeed such an extraordinary depth at these wavelengths was crucial to the principal 
aim of the HUDF09 program — namely, the search for z ~ 8-10 candidates appearing as Y105/J125 drop-outs (Bouwens 
et al. 2010a, 2011). Consequently, the contributions of both background sky noise and shot noise in the ERS are 
~a/14 times higher than in the HUDF09. The final image noise in our master frames is, moreover, weakly correlated 
between adjacent pixels as a result of the reduction process. Hence, rather than simply adding an independent noise 
component to each pixel to simulate the degradation in sky noise we attempt to at least partially approximate the 
impact of pixel-to-pixel correlation by degrading our UDF images instead through the stacking of blank regions of real 
background sky. For the addition of shot noise, however, we simply suppose that all flux detected above the 3<r sky 
noise level in the original UDF images corresponds to real object flux, and we directly resample this flux (with Poisson 
uncertainty) in the R statistical package at the ERS rate. 

In Figure ^ we compare the original (i?i6o) HUDF09 postage stamp images of three archetypal CLASS 1 (disk, 



or bulge/disk) systems from Fig. |8J and three archetypal CLASS 3 (irregular) systems from Fig. 10, against our new 
versions of these images degradea to the noise level of the ERS. It is clear from inspection of this Figure that the 
impact of the additional noise is to reduce the relative prominence of certain faint and/or fine-structure features in these 
z ~ 2 galaxies, but that the broad morphological properties of each system (e.g., the presence/absence of multiple 
nuclei, a central bulge, an outer disk, or a warped disk) remain readily discernable, and thus each galaxy remains 
clearly recognisable as a member of its assigned class. Of course, studies attempting to achieve finer morphological 
classifications than those conducted here, or else attempting to perform quantitative 'clump-science' (e.g., Forster- 
Schreiber et al. 2011), may well prove more sensitive to the imaging depth than our relatively crude analysis of Section 
[7| Further, one should not attempt to infer from the results of this Appendix that certain quantitative morphological 
analysis procedures (e.g., bulge-disk decomposition) will also prove insensitive to the difference in the HUDF09 and 
ERS imaging depths — such a statement clearly requires a more detailed investigation tailored to each morphological 
analysis tool under study. 

As a further robustness test of the morphological evolution described in Section[7|we have investigated the possibility 
that the decreasing fraction of massive spheroids with increasing redshift (past z ~ 2) is an artifact of the impact of 
cosmological surface brightness dimming (and band-pass shifting) on our observations. In particular, we have verified 
that our visual classifications of all massive elliptical galaxies in our 1.5 < z < 2.15 sub-sample remain consistent after 
artificial redshifting of their i^i6o images to z = 3. The results of this redshifting for six archetypal z ~ 1.5 ellipticals 
from our sample are displayed in Fig. [13] The extremely centrally-concentrated light profiles of these systems ensure 
that even under substantial cosmological surface brightness dimming they remain quite distinct in appearance from 
the remaining three classes. 
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Fig. 13. — Confirmation of the robustness of our visual morphological classifications of 1.5 < z < 3.5 galaxies in the ERS rest-frame optical 
(if 160 ) imaging of the GOODS-South field with WFC3 despite the greater contribution from background sky and shot noise relative to the 
deeper HUDF09 data. Specifically, the original ( if i6o) HUDF09 postage sta mp images of three archetypal CLASS 1 (disk, or bulge/disk) 
systems from Fig. [8] and three archetypal CLASS 3 (irregular) systems from Fig. |10| are presented in the top row of this Figure for comparison 
against versions of these images degraded to the noise level of the ERS in the bottom row. Although the faint morphological features of 
these galaxies appear less distinct after degradation they are still clearly recognisable as members of their assigned classes. Hence, we argue 
that the depth of the ERS imaging is indeed sufficient for confident morphological classification at the (crude) level of detail we employ in 
the analysis of Section [7] 




Fig. 14. — Confirmation of the consistency of visual morphological classifications of ellipticals at z ~ 1.5 and z ~ 3 in the ERS rest-frame 
optical (if 160 ) imaging of the GOODS-South field with WFC3. Specifically, the original (if i6o) ERS postage stamp images of six archetypal 
CLASS 2 spheroids at z ~ 1.5 are presented in the top row of this Figure for comparison against versions of these images artificially redshifted 
to z ~ 3 in the bottom row. The extremely centrally-concentrated light profiles of these giant ellipticals ensure that even under substantial 
cosmological surface brightness dimming they remain morphologically recognisable and quite distinct from the remaining three classes. 
Note that in the cases of ERS2182 and ERS 1238 we artificially redshift also their nearby companions, which our photometric (and/or 
spectroscopic) redshifts indicate to be at very nearly identical epochs. 



